The dark side of gravity: Modified theories of gravity Francisco S. N. LobcB Centro de Astronomia e Astrofisica da Universidade de Lisboa, Campo Grande, Ed. C8 1749-016 Lisboa, Portugal and Institute of Gravitation & Cosmology, University of Portsmouth, Portsmouth P01 2EG, UK (Dated: July 10, 2008) Modern astrophysical and cosmological models are faced with two severe theoretical difficulties, that can be summarized as the dark energy and the dark matter problems. Relative to the former, it has been stated that cosmology has entered a 'golden age', in which high-precision observational data have confirmed with startling evidence that the Universe is undergoing a phase of accelerated expansion. Several candidates, responsible for this expansion, have been proposed in the literature, in particular, dark energy models and modified gravity, amongst others. One is liable to ask: What is the so-called 'dark energy' that is driving the acceleration of the universe? Is it a vacuum energy or a dynamical field ("quintessence")? Or is the acceleration due to infra-red modifications of Einstein's theory of General Relativity? In the context of dark matter, two observations, namely, the behavior of the galactic rotation curves and the mass discrepancy in galactic clusters, suggest the existence of a (non or weakly interacting) form of dark matter at galactic and extra-galactic scales. It has also been proposed that modified gravity can explain the galactic dynamics without the need of introducing dark matter. We briefly review some of the modified theories of gravity that address these two intriguing and exciting problems facing modern physics. I. INTRODUCTION CT. u . bJO. > o r^ Cosmology is said to be thriving in a golden age, where a central theme is the perplexing fact that the Universe is undergoing an accelerating expansion [l|. The latter, one of the most important and challenging current problems in cosmology, represents a new imbalance in the governing gravitational equations. Historically, physics has addressed such imbalances by either identifying sources that were previously unaccounted for, or by altering the governing equations. The cause of this acceleration still remains an open and tantalizing question. The standard model of cosmology has favored the first route to addressing the imbalance, namely, a missing stress- y/S ■ energy component. In particular, the dark energy models arc fundamental candidates responsible for the cosmic expansion (see Refs. 2] for a review and references therein). A simple way to parameterize the dark energy is by an equation of state of the form ui = p/p, where p is the spatially homogeneous pressure and p the energy density of the dark energy. A value of uj < —1/3 is required for cosmic expansion, as dictated by the Friedmann equation a J a — ~AttG{p + p/3), and to = — 1 corresponds to a cosmological constant. A possibility that has been widely ^^ ■ explored, is that of quintessence, a cosmic scalar field cj) that has not yet reached the minimum of its potential V{4>) [3|]. A common example is the energy of a slowly evolving scalar field with positive potential energy, similar to the inflaton field used to describe the inflationary phase of the Universe. In quintessence models the parameter range is — 1 < lo < —1/3, and the dark energy decreases with a scale factor a(t) as pq ex a _3< - 1+w ^ J4j. A specific exotic form of dark energy denoted phantom energy, with w < — 1, has also been proposed [5(, and possesses peculiar properties, such as the violation of the energy conditions and an infinitely increasing energy density. However, recent fits to supernovae, cosmic microwave background radiation (CMBR) and weak gravitational lensing data indicate that an evolving equation_of state crossing the phantom divide, is mildly favored, and several models have been proposed in the literature [6, 7]. In particular, models considering a redshift dependent equation of state, possibly provide better fits to the most recent and reliable SN la supernovae Gold dataset. It is also interesting to test specific models that are motivated by particle physics against the SN data, rather than trying to fit the phenomenological equations of state. In a cosmological setting, it has also been shown that the transition into the phantom regime, for a single field is probably physically implausible [6j, so that a mixture of various interacting non-ideal fluids is necessary. If confirmed in the future, this behavior has important implications for theoretical models of dark energy. For instance, this implies that dark energy is dynamical and excludes the cosmological constant and the models with a constant parameter as possible candidates for dark energy. An alternative model to dark energy is that of the generalized Chaplygin gas (GCG), based on a negative pressure fluid, which is 'Electronic address: flobo@cosmo.fis.fc.ul.pt inversely proportional to the energy density [8||, i.e., p c h = —A/p" h , where A and a are positive constants. An attractive feature of this model, is that at early times, the energy density behaves as matter, p^ ~ a~ 3 , where a is the scale factor, and as a cosmological constant at a later stage, p c h = const. This dual behavior is responsible for the interpretation that the GCG model is a candidate of a unified model of dark matter and dark energy [9], and probably contains some of the key ingredients in the dynamics of the Universe for early and late times. All of these models present an extremely fascinating aspect for future experiments focussing on supernovae, CMBR and weak gravitational lensing and for future theoretical research. One may also explore the alternative viewpoint, namely, through a modified gravity approach. A very promis- ing way to explain these major problems is to assume that at large scales Einstein's theory of General Relativity breaks down, and a more general action describes the gravitational field. The Einstein field equation of General Relativity was first derived from an action principle by Hilbert, by adopting a linear function of the scalar curva- ture, R, in the gravitational Lagrangian density. However, there are no a priori reasons to restrict the gravitational Lagrangian to this form, and indeed several generalizations of the Einstein-Hilbert Lagrangian have been proposed, including "quadratic Lagrangians" . involving second order curvature invariants such as i? 2 , R^R^, R a f3 l j,vR al3 ^ ', e al3 ^ v Rap-ysKly, Capfj.uC ' 11 ^ 1 ' , etc [1J|. The physical motivations for these modifications of gravity were related to the possibility of a more realistic representation of the gravitational fields near curvature singularities and to create some first order approximation for the quantum theory of gravitational fields. In this context, a more general modifi- cation of the Einstein-Hilbert gravitational Lagrangian density involving an arbitrary function of the scalar invariant, f(R), was considered in [ll| . and further developed in [12, [Da, LLJ] • Recently, a renaissance of f(R) modified theories of gravity has been verified in an attempt to explain the late-time accelerated expansion of the Universe. Earlier interest in f(R) theories was motivated by inflationary scenarios as for instance, in the Starobinsky model, where f(R) = R — A + aR 2 was consider ed fl5l| . In particular, it was shown that cosmic acceleration can be indeed explained with the context of f(R) gravity [la ], and the conditions of viable cosmological models have also been derived [171 ]. In the context of the Solar System regime, severe weak field constraints seem to rule out most of the models proposed so far [H, [l9| , although viable models do exist [20] • In the context of dark matter, the possibility that the galactic dynamics of massive test particles may be understood without the need for dark matter was also considered in the framework of f(R) gravity models [2l|, [22|, [23, [2J, |2a |. and connections with MOND and the pioneer anomaly further explored by considering an explicit coupling of an arbitrary function of R with the matter Lagrangian density [261 |27| . The issue of dark matter is a long outstanding problem in modern astrophysics. Two observational aspects, namely, the behavior of the galactic rotation curves and the mass discrepancy in clusters of galaxies led to the necessity of considering the existence of dark matter at a galactic and extra-galactic scale. The rotation curves of spiral galaxies show that the rotational velocities increase from the center of the galaxy and then attain an approximately constant value, Vt goo ~ 200 — 300 km/s, within a distance r from the center of the galaxy [28[ . In these regions the mass increases linearly with the radius, even where very little luminous matter can be detected. Relatively to the mass discrepancy in clusters of galaxies, the total mass of a cluster can be estimated in two ways. First, by taking into account the motions of its member galaxies, the virial theorem provides an estimate, My. Second, the total baryonic mass M may be estimated by considering the total sum of each individual member's mass. The mass discrepancy arises as one generally verifies that My is considerably greater than M, with typical values of My jM ~ 20 — 30 [28[. This is usually explained by postulating the existence of a dark matter, assumed to be a cold pressure-less medium distributed in a spherical halo around the galaxies. Still in the context of modified gravity, an interesting possibility is the existence of extra dimensions. It is widely believed that string theory is moving towards a viable quantum gravity theory, and one of the key predictions of string theory is precisely the existence of extra spatial dimensions. In the brane-world scenario, motivated by recent developments in string theory, the observed 3-dimensional universe is embedded in a higher-dimensional spacetime [29j . Most brane-world models, including those of the Randall- Sundrum type [30j, produce ultra-violet modifications to General Relativity, with extra-dimensional gravity dominating at high energies. However it is also possible for extra- dimensional gravity to dominate at low energies, leading to infra-red modifications of General Relativity. New features emerge in the brane scenario that may be more successful in providing a covariant infra-red modification of General Relativity, where it is possible for extra-dimensional gravity to dominate at low energies. The Dvali-Gabadadze-Porrati (DGP) models [3l| achieve this via a brane induced gravity effect. The generalization of the DGP models to cosmology lead to late-accelerating cosmologies [32| , even in the absence of a dark energy field 33]. This exciting feature of "self acceleration" may help towards a new resolution to the dark energy problem, although this model deserves further investigation as a viable cosmological model [34[. While the DGP braneworld offers an alternative explanation to the standard cosmological model, for the expansion history of the universe, it offers a paradigm for nature fundamentally distinct from dark energy models of cosmic acceleration, even those that perfectly mimic the same expansion history. It is also fundamental to understand how one may differentiate this modified theory of gravity from dark energy models. The DGP braneworld theory also alters the gravitational interaction itself, yielding unexpected phcnomcnological extensions beyond the expansion history. Tests from the solar system, large scale structure, lensing all offer a window into understanding the perplexing nature of the cosmic acceleration and, perhaps, of gravity itself [34]]. The structure formation (35[ and the inclusion of inflation are also important requirements of DGP gravity, if it is to be a realistic alternative to the standard cosmological model. Generalizations of the DGP model with the inclusion of a Gauss-Bonnet (GB) term have also been explored [33], the global structure of the DGP cosmologies have been analyzed [36J , and research into other brane- world approaches to dark energy, such as the supersymmetric large extra dimensions (SLED) model have been undertaken [371 ]. In this work, we review several modified theories of gravity exploring some of their interesting properties and characteristics, in particular, focussing mainly on the late-time cosmic acceleration. We refer the reader to excellent reviews, for instance, on dark energy and modified gravity in Ref. Q, f(R) gravity in Ref. [38|, an introduction to modified gravity as an alternative for dark energy in Ref. [39(, and the late-time acceleration in braneworlds in Ref. [4fJ. This paper is outlined in the following manner. In Section [HI we review f(R) modified theories of gravity, for instance, focussing on the scalar-tensor representation of f(R) and on the late-time acceleration; we further analyze an interesting extension to f{R) gravity by considering an R— matter coupling, which yields some intriguing properties; we also discuss Gauss-Bonnet gravity and the late-time acceleration, and briefly consider modified Gauss-Bonnet gravity; the DGP brane-world model is also briefly reviewed. In Section [TTl] we consider the possibility of dark matter being a geometric effect in f(R) gravity, focussing on the f(R) generalized virial theorem, and its applications to galactic cluster observations. Finally in Section [IV! we conclude with a summary and discussions. II. MODIFIED THEORIES OF GRAVITY: LATE-TIME COSMIC ACCELERATION A. f(R) modified theories of gravity A promising avenue that has been extensively investigated recently are the f(R) modified theories of gravity, where the standard Einstein- Hilbert action is replaced by an arbitrary function of the Ricci scalar R J4ll |. In this work, we use the metric formalism, which consists in varying the action with respect to q^", although other alternative approaches have been considered in the literature, namely, the Palatini formalism [42|, [43|, where the metric and the connections are treated as separate variables; and the metric-affine formalism, where the matter part of the action now depends and is varied with respect to the connection [431 ] . 1. Action and field equations The action for the f(R) modified theories of gravity is given by S=^Jd i xV^f(R) + S M (9^^), (1) where k = 8ttG. Sm(<7' w ') ip) is the matter action, defined as Sm = / d^x^J—g C m (g^ u ,tlj), where C m is the matter Lagrangian density, in which matter is minimally coupled to the metric g^ v and ip collectively denotes the matter fields. Using the metric approach, by varying the action with respect to g 1 *", provides the following field equation FR»» - \f g„u - V M V„F + g^DF = k t£?> , (2) where F = df/dR. The matter stress-energy tensor, T^™, is defined as "" V~9 S(g^) ■ l > Now, considering the contraction of Eq. ^\j, provides the following relationship FR-2f + 3DF= kT, (4) which shows that the Ricci scalar is a fully dynamical degree of freedom. Note that the field equation, Eq. ([5]), may be written as 1 cff (j>l/ = Rp,U — nRd^U = K ^ fill J (5) where the effective stress energy tensor is given by T°^ = TJfJ + Tfi? ■ The components TJu? and the curvature stress-energy tensor, T^V , are defined as fjZ? = TW/F, (6) T$ = -^ V IM V„F-~g tiV (RF + nF + K T) (7) respectiveiy. It is also interesting to consider the conservation law for the above curvature stress-energy tensor. Taking into account the Bianchi identities, V M G M „ = 0, and the diffeomorphism invariance of the matter part of the action, which yields V M T^™ = 0, we verify that the effective Einstein field equation provides the following conservation law V^t( c ) = rHycp (r\ 2. Scalar-tensor representation for f(R) gravity f(R) gravity may be written as a scalar-tensor theory, by introducing a Legendre transformation {R, /} — > {(f>, V} defined as cj> = F(R), V(<f>) = R(<P)F-f(R(d>)). (9) In this representation the field equations of f(R) gravity can be derived from a Brans-Dicke type action with parameter u> = 0, given by S=^~ f[<f>R- V{4>) + L m ] JS d 4 x. (10) 2k J The only requirement for the gravitational field equations to be expressed in the form of a Brans-Dicke theory is that F(R) be invertible, that is, R(F) exists p9{ . Thus, the field equations of f(R) gravity can be reformulated as 1 G Rixv - -^g^R = 8tt—t^ + 0^ v , (n) where V = ~\v [4>) 9nu + -r (V M V, - 9ltv U) cj>. (12) Using these variables the trace equation, i.e., Eq. ([5]), takes the form 3D<£ + 2V{6) - (f>-rr = kT. (13) dep The modification of the standard Einstein-Hilbcrt action leads to the appearance of an effective gravitational constant G c g = G/(j> in the field equations. Note the presence of a new effective source term for the gravitational field, given by the tensor 6^. 3. R— matter couplings in f(R) gravity Recently, in the context of f(R) theories of modified gravity, it was shown that a function of R— matter coupling induces a non-vanishing covariant derivative of the stress-energy, V ^F" ^ 0. This potentially leads to a deviation from geodesic motion, and consequently the appearance of an extra force [261 ]. Implications, for instance, for stellar equilibrium have been studied in Ref. [27[ ■ The equivalence with scalar-tensor theories with two scalar fields has been considered in Ref. [4J], and a viability stability criterion was also analyzed in Ref. [45|. It is interesting to note that nonlinear couplings of matter with gravity were analyzed in the context of the accelerated expansion of the Universe [46J , and in the study of the cosmological constant problem [47 1 . The action for R— matter couplings, in f(R) modified theories of gravity [261 ] . takes the following form 1 S = -f 1 (R) + [l + Xf 2 (R)]C 7 -g d x (14) where fi(R) (with i — 1,2) are arbitrary functions of the curvature scalar R. For notational simplicity we consider k — 1 throughout this subsection. Varying the action with respect to the metric g^ yields the field equations, given by FiR„ v ~ \h9y.v ~ V M V, Fx + g^UFx = -2XF 2 C m R^ + 2A(V„V„ - g^U)C m F 2 + (1 + Xf 2 )T^ , (15) where we have denoted Fi(R) = fl(R), and the prime represents the derivative with respect to the scalar curvature. Now, taking into account the generalized Bianchi identities [26|, [48| , one deduces the following corrected conservation equation V M T (m) = XF2 ^ l + A/ a a r - T( m ) VR, (16) where the coupling between the matter and the higher derivative curvature terms describes an exchange of energy and momentum between both. In the following, consider the equation of state for a perfect fluid T^ = {p + p)U,U„ + pg^u (17) where p is the energy density and p, the pressure, respectively. The four- velocity, U^, satisfies the conditions U^U^ = -1 and U^U^u = 0. Introducing the projection operator /i M „ = g M „ + U^U V , gives rise to non-geodesic motion governed by the following equation of motion for a fluid element: dU 1 * jds + T^plf^U 13 — / M , where the extra force, / M , is given by f" P XF 2 1 + A/j {Z m -p)VuR + V v p h^ (18) In a recent paper 49], the authors argued that a "natural choice" for the matter Lagrangian density for perfect fluids is C m = p, based on Refs. [5Ctl5l|. where p is the pressure. This choice has a particularly interesting application in the analysis of the R— matter coupling for perfect fluids, which implies in the vanishing of the extra force [261 ]. However, it is important to point out that despite the fact that C m — p does indeed reproduce the perfect fluid equation of state, it is not unique [52J ]. Other choices include, for instance, C m = — p [5X|. l53||. where p is the energy density, or C m — ~na, were n is the particle number density, and a is the physical free energy defined as a = p/n — Ts, with T being the temperature and s the entropy per particle (see Ref. [5ll.l52j] for details). Hence, it is clear that no immediate conclusion may be extracted regarding the additional force imposed by the non-minimum coupling of curvature to matter, given the different available choices for the Lagrangian density. One may conjecture that there is a deeper principle or symmetry that provides a unique Lagrangian density for a perfect fluid [53]. This has not been given due attention in the literature, as arbitrary gravitational field equations depending on the matter Lagrangian have not always been the object of close analysis. See Ref. [52] for more details. 4. Late-time cosmic acceleration In this subsection, we show that f(R) gravity may lead to an effective dark energy, without the need to introduce a negative pressure ideal fluid. Consider the FLRW metric ds l -dt 2 + a 2 (t) dr 2 1 — kr 2 + r 2 (d6 2 + sin 2 9d(f> 2 ) (19) Taking into account the perfect fluid description for matter given by Eq. (|17p . we verify that the gravitational field equation, Eq. ([5]), provides the generalised Friedmann equations in the following form [54ll55|: 1 2F(R) {a ^{llf(K)-mR)}-Z^)RF>(R)} = Ip, 1 -RF'{R) + RF'(R) + R 2 F"(R)--[f(R)-RF{R:)}} = --(p + 3p) 6 (20) (21) These modified Friedmann field equations may be rewritten in a more familiar form, as a \ k a J 6 -) = ~(ptot+3ptot), (23) a J b where p tot — p + pt c \ and p tot = p + Pm, and the curvature stress-energy components, pt c \ and pt c )i are defined as P (c) = J^^im-RFm-S^RF'iR)} , (24) P {c) = ^y J2 Q RF'(R) + RF'(R) + R 2 F'\R) - \ [f(R) RF(R)]j , (25) respectively. The late-time cosmic acceleration is achieved if the condition p to t + 3ptot < is obeyed, which follows fromEq. pg|) . For simplicity, consider the absence of matter, p = p = 0. Now, taking into account the equation of state w e ff — P(c)/p(c)i with f(R) oc R n and a generic power law a(t) = ao(£/£o) a [54j ]. the parameters w g and a are given by 6n 2 - In - 1 -2n 2 + 3n - 1 Wpff = ^ , a = , (26 6n 2 - 9n + 3 ' n - 2 ' V ; respectively, for n 7^ 1. Note that a suitable choice of n can lead to the desired value of u> e g < —1/3, achieving the late-time cosmic acceleration. Another example is a model of the form f(R) = R — /i 2 ( n+1 ) / ' R n [lq . Choosing once again a generic power law for the scale factor, the parameter can be written as „ C „ = - 1 + _&±H> , 27) 6(2n + l)(n + 1) and once again a desired value of w e ff < — 1/3 may be attained, by appropriately choosing the value of the parameter n. Note that as n — * 00 the spacetime is approximately de Sitter. Other forms of f(R) have also been considered in the literature, for instance those involving logarithmic terms, such as /(i?) = R + a\n{R/n 2 )+0R m or f(R) = R + jR- n [\n(R/p 2 )] m [H, HI] . These models also yield acceptable values for the effective equation of state parameter, resulting in the late-time cosmic acceleration. B. Gauss-Bonnet gravity and cosmic acceleration 1. Gauss-Bonnet gravity In considering alternative higher-order gravity theories, one is liable to be motivated in pursuing models consistent and inspired by several candidates of a fundamental theory of quantum gravity. In this context, it may be possible that unusual gravity-matter couplings predicted by string/M-theory may become important at the recent low-curvature Universe. For instance, one may couple a scalar field not only with the curvature scalar, as in scalar-tensor theories, but also with higher order curvature invariants. Indeed, motivations from string/M-theory predict that scalar field couplings with the Gauss-Bonnet invariant Q are important in the appearance of non-singular early time cosmologies. It is also possible to apply these motivations to the late-time Universe in an effective Gauss-Bonnet dark energy model m- Consider the action of Gauss-Bonnet gravity given by j4 / ^ - \d^d^ - V(4>) + f{<t>)G ■SuisT, ■>!>), (28) where A = +1 is defined for a canonical scalar field, and A = — 1 for a phantom field, respectively. The Gauss-Bonnet invariant in given by Q = R 2 — AR^R^ + R^apR^ 01 ^ ' ■ It is also important to note that in the matter action, matter is minimally coupled to the metric and not to the scalar field, making Gauss-Bonnet gravity a metric theory. Thus, using the diffeomorphism invariance of Sm (q^ j V 1 ) yields the covariant conservation of the stress-energy tensor, V' J TJ'™ ) = 0. Varying the action with respect to 0, provides the equation of motion for the scalar field, given by \V 2 0-V'( ( f>)+f'(0)g = Q. (29) The gravitational field equations are given by varying the action with respect to the metric g pv ', and provides the following relationship [55|, [57| -G^ - \gTf{<t>)Q + 2f(cj>)RR^ - 2WV v [f(cf>)R} + 2g pv V 2 [f{4>)R] K 2 ~8f(cf>)R%R»? + 4V p VTO)iH + iVpVmW"] - 4V 2 [/(0)i?H -Ag^V p V a [fWR P l + 2f{$)R ppaT R% aT - 4V P V CT [/(<£) R^ v ] = T pv + T$ v , (30) where Tj? v is given by Id^Q^-lg^d^dpA-l, TT = A ( tt^ d u 4> - -y v d p ^ d p <p ) - -g^V(4>) . (31) One may use the following identities obtained from the Bianchi identity [55l . 157] : VRpr^v = V„R vr ~ V„iV , Vi^ M = - V M i? , V /i V„i2'"' = -OR , V p V a R ppl/a = V 2 R^ - ^VViJ + R ppva R pa - R p p R vp , VpV (M^)P = v^V^-R - RWRpg + R p p R vp , in the Gauss-Bonnet gravitational field equation, which may then be formally simplified to K 2 +2f{cj>)R ppaT R l/ pcTT - 4f(0)R pp ™R pa = T pv + T$ v + T? v , (32) -G pv - ^g pv I{cp)Q + 2f(4>)RR^ + Af{<t>)R» p R v P where T pv is defined as T pv = 2[V p V v J{(t>)]R-2g pu [V 2 f{4>)]R-A[V p V p f{4>)]R vp - A[V p V u f{4>)]R pp +4{V 2 f(cb)}R p » + 4 5 ^[V p V ff /(0)]i?" CT - 4[V p V a f(4>)}R pp ™ . (33) Using the FLRW metric, the modified Friedmann equations for Gauss-Bonnet gravity reduce to the following relationships a \ k , k a J 6 3 d (35) A0 2 - 1/(0) + l2H^f{cj>) + 12- (ff 2 /) and the equation of motion for the scalar field is given by A U + ZH^j + V'{4>) - 2Af'{<t>)H 2 [H + H 2 ) = , (36) where the expression for the Gauss-Bonnet invariant, Q = 24H 2 (H + H 2 ), is used. Note that in the absence of matter, the Gauss-Bonnet gravitational field equations may be written as p GB - -H 2 , pgb = -- (SH 2 + 2H) , (37) where the Gauss-Bonnet curvature stress-energy tensor components are defined as Pgb = \tf + V{ct>) - 240/'(0)ff 3 , (38) Pgb = ^ 2 -y(^)+8^(i/ 2 /)+16ff 3 0/'(0). (39) These relationships are particularly interesting as one may now define an effective equation of state given by PGB 1 2H Pgb (40) Consider the choices of an exponential scalar-GB coupling and exponential scalar potential given by [57| V(cf>) = Voe-W* , f{4>) = foe 2 ^ . (41) The scale factor and the scalar field are respectively chosen in the following form a(t) a t h °, for ho > (quintessence), ao(t s — t) h ° , for ho < (phantom), and m .In for fto > 0, ^o In ( %-* ) , for /i < 0, where ti is an arbitrary constant [57|. Using the above choices, and considering the absence of matter, then taking into account Eqs. finally reorganizing the terms leads to the following equations Vot\ 48/ feg 1 k(1 + M 6 k(1 + /io) 3ftg(l - fto) A0qk(1 — 5/iq) /jo A0q« (42) (43) and ([3"6")l , and (44) (45) The effective equation of state parameter, Eq. (|4"0")) . takes the following form Weff = -1 - — . 3/l (46) Note that if /iq < then w e ff < — 1, reflecting an effective phantom regime; and if ho > then u> e & > — 1, which reflects an effective quintessence regime. However, it has been shown that the case of ho < is always stable, while the case of a non-phantom ho > cosmology is always unstable [57} • An interesting example is the case of Vo = 0, which from Eq. (|4~4"|) imposes the following condition [57[: </>g = — 6/Iq(1 — h )/[An(l — 5h )]. In order for cf> to be real, in the case of the a canonical scalar, A = 1, one finds that 1/5 < ho < 1; for a phantom field, A = — 1, then /io < 1/5 or /io > 1. To achieve an effective equation of state parameter value that mimics dark energy is to consider, for instance, ho = —80/3, which imposes the w c g = —1.025. We refer the reader to Ref. [57| for more examples and details. It is also interesting to note that observational constraints [58| , such as CMBR, galaxy distribution, large scale structure and supernovae seem to favour the Gauss- Bonnet coupling. 2. Modified Gauss-Bonnet gravity and the late-time acceleration An interesting alternative gravitational theory is modified Gauss-Bonnet gravity, which is given by the following action: S = -9 £ + w SmGT,^). (47) This theory has been extensively analyzed in the literature [39|, [59|, |6CJ , and rather than review all of its intricate details here, we note that it is a subset of Gauss-Bonnet gravity given by the action (j2"5|) . To see this, we follow closely the approach outlined in Ref. [5a, [60]. By introducing two auxiliary scalar fields A and B, the gravitational part of the action ([47]) . may be rewritten as S = j4 / £- + B(g-A) + f(A) (48) Now varying with respect to B, one obtains A = Q, so that the action (f48|) is recovered. Varying with respect to A, one obtains B = f'(A), and substituting in (|4"5)l leads to S = Id Xy/—g With the following definitions <p = A , and one finally ends up with the following action R R 2k + f'(A)(g-A) + f(A) V{<j>) = Af (A) - /(A) d xy—g 2 k V(^) + /($G s M (ar,ii>) (49) (50) (51) which is simply the action for Gauss-Bonnet gravity given by Eq. (j2"8|) with the absence of the kinetic term. Thus, modified Gauss-Bonnet theory, given by (|4"T|) is dynamically equivalent to Gauss-Bonnet gravity with A = |55l . l6(j | . We refer the reader to Refs. [39j, [59j for more details. C. DGP brane gravity and self-acceleration One of the key predictions of string theory is the existence of extra spatial dimensions. In the brane-world scenario, motivated by recent developments in string theory, the observed 3-dimensional universe is embedded in a higher- dimensional spacetime [29]. One of the simplest covariant models is the Dvali-Gabadadze-Porrati (DGP) braneworld model, in which gravity leaks off the AD Minkowski brane into the 51? bulk at large scales. The generalization of the DGP models to cosmology lead to late-accelerating cosmologies [32|], even in the absence of a dark energy field [33j |. While the DGP braneworld offers an alternative explanation to the standard cosmological model, for the expansion history of the universe, it offers a paradigm for nature fundamentally distinct from dark energy models of cosmic acceleration, even those that perfectly mimic the same expansion history. The 5D action describing the DGP model is given by S 1 2k 5 cPxy/=g w R (5), 1 2^4 d A x, Mr- d X\f— 7 C r (52) where k$ = 8nG5 and K4 = %-kG^. The first term in the action is the Einstein- Hilbert action in five dimensions for a five-dimensional bulk metric, gAB, with a five-dimensional Ricci scalar ^R; the second term is the induced Einstein-Hilbert term on the brane, with a four-dimensional induced metric 7 on the brane; and C m represents the matter Lagrangian density confined to the brane. The transition from AD to hD behavior is governed by a cross-over scale, r c , given by ^5 2K4 (53) Gravity manifests itself as a 4-dimensional theory for characteristic scales much smaller than r c ; for large distances compared to r c , one verifies a leakage of gravity into the bulk, consequently making the higher dimensional effects important. Thus, the leakage of gravity at late times initiates acceleration, due to the weakening of gravity on the brane. For weak fields, the gravitational potential behaves as $ for for r < r c , r > r c . (54) by Taking into account the FRW metric, and considering a flat geometry, the modified Friedmann equation is given H 2 -H lirG -p. (55) where e — ±1, and the energy density satisfies the standard conservation equation, i.e., p + 3H(p + p) = 0. For scales H^ 1 <C r c , the second term is negligible, and Eq. ([55)1 reduces to the general relativistic Friedmann equation, 10 i.e., H 2 = 8irGp/3. The second term becomes significant for scales comparable to the cross-over scale, H~ 1 > r c . Self- acceleration occurs for the branch e = +1, and the modified Friedmann equation shows that at late times in a CDM universe characterized by a scale factor p ex a~ 3 , the universe approaches a de Sitter solution H^H oc = - (56) r c Thus, one may achieve late-time acceleration if the Hq is of the order of ffoo. Note that the late-time acceleration in the DGP model is not due to the presence of a negative pressure, but simply due to the weakening of gravity on the brane as a consequence of gravity leakage at late times. Although the weak-field gravitational DGP behaves as AD on scales smaller than r c , linearized DGP gravity is not described by General Relativity [4CJ, [61j . It is also interesting to note that despite the fact that the expansion history of the DGP model and General Relativity are the same, the structure formation in both are essentially different [62]. Combining these features provides the possibility of distinguishing the DGP model from dark energy models in General Relativity. We should also emphasize that while the DGP model is not ruled out by current observations, the ACDM model fits the data comfortably [63|] . Another interesting aspect of this model is that the self-accelerating branch in the DGP model contains a ghost at the linearized level [6l|, [64]]. The presence of the ghost implies a negative sign for the kinetic term, resulting in negative energy densities, consequently leading to the instability of the spacetime. However, in a recent paper it was claimed that a higher codimension generalization of the DGP scenario is free of ghost instabilities [63 , and further work along these lines is currently underway. We refer the reader to Ref. (401. l6lj. and references therein, for more details on the DGP model. III. DARK MATTER AS A GEOMETRIC EFFECT OF MODIFIED GRAVITY The issue of dark matter is a long outstanding problem in modern astrophysics. Two observational aspects, namely, the behavior of the galactic rotation curves and the mass discrepancy in clusters of galaxies led to the necessity of considering the existence of dark matter at a galactic and extra-galactic scale. The galactic rotation curves of spiral galaxies 28] are probably the most striking evidences for the possible failure of Newtonian gravity and of General Relativity on galactic and intergalactic scales. In these galaxies, neutral hydrogen clouds are observed at large distances from the center, much beyond the extent of the luminous matter. As these clouds are moving in circular orbits with nearly constant tangential velocity v tg , such orbits are maintained by the balance between the centrifugal acceleration v 2 /r and the gravitational attraction GM(r)/r 2 of the total mass M(r) contained within the radius r. This yields an expression for the galactic mass profile of the form M(r) = rv 2 /G, with the mass increasing linearly with r, even at large distances, where very little luminous matter has been detected [281 ] . This peculiar behavior of the rotation curves is usually explained by postulating the existence of dark matter, assumed to be a cold and pressureless medium, distributed in a spherical halo around the galaxies. There are many possible candidates for dark matter, the most popular ones being the weakly interacting massive particles (WIMP) [69 ]. One cannot also a priori exclude the possibility that Einstein's (and Newton's) theory of gravity breaks down at galactic scales. In this context, several theoretical models, based on a modification of Newton's law or of General Relativity, have been proposed to explain the behavior of the galactic rotation curves [671 ] . A promising avenue that has been extensively investigated recently are the f(R) modified theories of gravity. In this context, early work in explaining dark matter using f(R) gravity using models of the form f{R) ex R n found large values for n [2l|, [22j, [23j • In these papers, a power law modified Newtonian potential of the form <f>(r) = — ^?[1 + {r/r c )P] was considered, to describe the observed behavior of the galactic rotation curves, where m is the mass of the particle, r c a constant and the coefficient (3 depends on the 'slope' parameter n in the modified action. Using this modified Newtonian potential, it was found that the best fit to 15 low luminosity rotation curves in R n gravity is obtained for n — 3.5 |2l| (somewhat lower values, in particular, n — 2.2, were obtained in [22, UM)- These results seem to suggest that a strong modification of standard general relativity is required to explain the observed behavior of the galactic rotation curves. Note that these large values of n are in gross violation with the Solar System tests. However, recently is was found that only slight deviations from General Relativity are needed, i.e., n = 1 + e with e <C 1 [2J, |2J|. This discrepancy of values can be traced back to the correction term of the modified Newtonian potential. It was shown in Ref. [24| that the correct modified term to the Newtonian potential in the "dark matter" dominated region, where the rotation curves are strictly flat, must have a logarithmic dependence on the radial coordinate r, of the form $Ar(r) = — ^^ + v 2 g ln(r/ro), where ro is an arbitrary constant of integration (we refer the reader to Ref. 24] for details). These differences in the Newtonian limit in the two models result in different values of the parameter n in the power-law modified action of the gravity. In the following sections, we shall analyze the 'dark matter' problem by considering a generalized version of the virial theorem in the framework of f(R) modified theories of gravity [25[. Recall that due to its generality and wide 11 range of applications, the virial theorem plays an important role in astrophysics. Assuming steady state, one of the important results which can be obtained with the use of the virial theorem is to deduce the mean density of astrophysical objects such as galaxies, clusters and super clusters, by observing the velocities of test particles rotating around them. Hence the virial theorem can be used to predict the total mass of the clusters of galaxies. The generalized virial theorem, in the context of f(R) gravity is obtained by using a method based on the collisionless Boltzmann equation [25| . The additional geometric terms present in the modified gravitational field equations provide an effective contribution to the gravitational energy, which at the galactic/extra-galactic level acts as an effective mass, playing the role of the 'dark matter'. The total virial mass of the galactic clusters is mainly determined by the effective mass associated to the new geometrical term, the geometrical mass. It is important to note that the latter term may account for the well-known virial theorem mass discrepancy in clusters of galaxies. A. The generalized virial theorem in f(R) gravity Consider an isolated and spherically symmetric cluster described by a static and spherically symmetric metric ds 2 = -e^dt 2 + e x ^dr 2 + r 2 {dd 2 + sin 2 ^ 2 ) . (57) The galaxies, treated as identical and collisionless point particles, are described by a distribution function fg, which obeys the general relativistic Boltzmann equation. In terms of the distribution function the stress-energy tensor can be written as T^ v = J f B mu^u^ du, where m is the mass of the particle (galaxy) [68j |. u^ is the four- velocity of the galaxy and du = du r duedu v / u t is the invariant volume element of the velocity space. Thus, the stress-energy tensor of the matter in a cluster of galaxies can be (r) represented in terms of an effective density p c g and of an effective anisotropic pressure, with radial p y c g and tangential p cS components, given by Pcti = P (u 2 ) , p^ = p (u 2 r ) , p[^ ] = p (uf) = p (u 2 ,) , (58) where, at each point, (w 2 ) is the average value of u 2 , etc, and p is the mass density [69J | . In what follows, we use this form of the stress-energy tensor, and for convenience take into account the scalar- tensor representation of f(R) gravity outlined in Section Hi A 21 As we are interested in astrophysical applications at the extra-galactic level, we may assume that the deviations from standard General Relativity (corresponding to the background value <j> = 1) are small. Therefore we may represent as = 1 4- eg'(R), where e is a small quantity, and g'(R) describes the modifications of the geometry due to the presence of the tensor 9^ v [19|], so that 1/0 ~ 1 — eg'(R). Now adding up the non-zero components of the gravitational field equation Eq. (fTTj) . and taking into account the above approximations (see Ref. [25j for details), one obtains the following relationship / a ii i i/'\'\ er x 1^— + — + -- — J ~ AnGp (u 2 ) + ^G P4n (59) where (u 2 ) = (u 2 ) + (u 2 ) + (w 2 ,) + (u 2 ), and the useful quantity p^ is defined as -ep(u 2 )g'(R) + ' 4ttG \v {(j))+ i(2V t V* + D) (60) b=l+eg'(R) which may be interpreted as the geometric energy density. It is convenient to introduce some approximations that apply to test particles in stable circular motion around galaxies, and to the galactic clusters. First of all, we assume that v and A are slowly varying (i.e. v' and A' small), so that in Eq. ([5"9"]l the quadratic terms can be neglected. Secondly, we assume that the galaxies have non-relativistic velocities, so that (u 2 ) « (u 2 ,) « (u§) <C (uf,) w 1. Thus, Eq. (|55j) becomes ±1 ( r ^ 2r 2 dr \ dr AnGp + 4irG P4> . (61) In order to derive the virial theorem for galaxy clusters, one uses the relativistic Boltzmann equation, which provides the following relationship (see Ref. [25| for details) 2K-- 4Trr 3 p^-dr = 0, with K = 2np [(uj) + (uj) + (u 2 )] r 2 dr . (62) 2 Jo " r Jo 12 K is the total kinetic energy of the galaxies, and the total mass of the system is given by M = /„ dM(r) = L 4wpr 2 dr. The main contribution to M is due to the baryonic mass of the intra-cluster gas and of the stars, but other particles, such as massive neutrinos, may also contribute significantly to M. Now, multiplying Eq. (|6ip by r 2 and integrating from to r we obtain GM{r) = -r 2 -^-~GM^(r), with Af (r) = 4tt / P4> {r')r n dr' . (63) The useful quantity M^ is denoted as the geometric mass of the cluster. By multiplying Eq. (f6"3"|) with dM(r), followed by an integration one deduces the relationship n = n^--j i^p-^dr, (64) with the following definitions n=-[ R ^n dM{r) , and n^^^MdMir), (65) o r where the quantity fl is the usual gravitational potential energy of the system. Finally, with the use of Eq. (I62|) . we obtain the generalization of the virial theorem, in f(R) modified theories of gravity, which takes the form IK + ft - ft = 0. (66) Note that the generalized virial theorem, given by Eq. (166|) , can be written in an alternative form if we introduce the radii Ry and R^ defined by Ry=M 2 / fMWdMM, and R+ = Mil C ^^ dM(r) , (67) respectively. We denote i?^ as the geometric radius of the cluster of galaxies. Thus, the quantities f2 and Q<p are finally given by O = —GM 2 /Ry and Cl^ = GMl/R^, respectively. The virial mass My is defined as 2K = ^^ . (68) Rv After substitution into the virial theorem, given by Eq. (|66[) , we obtain M + M 2 R 4> • [ > If My/M > 3, a condition which is valid for most of the observed galactic clusters, then Eq. (j6"5|) provides the virial mass in f(R) gravity, which can be approximated by Mv ^^L. (70) M R^ V ' B. Geometric mass and geometric radius from galactic cluster observations An interesting application of the generalized virial theorem can be inferred from the galaxy cluster observations. According to the modified f(R) gravity model, the total mass of the cluster consists of the sum of the baryonic mass (mainly the intra-cluster gas), and the geometric mass, so that Af tot (r) = in J„ (p g + p$) r 2 dr. Hence it follows that Mtot(f) satisfies the following mass continuity equation dM 2 (r) =^r 2 p g (r)+^r 2 PlP {r). (71) 13 Note that most of the baryonic mass in the clusters of galaxies is in the form of the intra-cluster gas. The gas mass density p g distribution can be fitted with the observational data by using the following expression for the radial baryonic mass (gas) distribution [70] / 2 \ -3/3/2 Pg(r) = Po ( 1 + ^J , (72) where r c is the core radius, and po and_5 are (cluster-dependent) constants. Using the Jeans equation [28J|, one may obtain the total mass distribution [7(J, [71( s so that taking into account the density profile of the gas given by Eq. (|72[) , the total mass profile inside the cluster is given by M tot (r)= 3kBPT J J* , (73) pm p G r£ + r* where ks is Boltzmann's constant, T g is the gas temperature, p rs 0.61 is the mean atomic weight of the particles in the cluster gas, and m p is the proton mass [701 ]. Using the mass continuity equation, Eq. ([TTj) . then Eqs. (172p and f73"|) provide the expression of the geometric density term inside the cluster, given by 3k B 0T g (r 2 +3r 2 c ) _ AnGpo 47rp0 (r) In the limit r > ^ we obtain for p^ the simple relation ^^ - AnGporfr 2 -^} ' and the geometric mass in the limit r » r c , may be approximated as ~3k B (3T g inGpor^r 2 - 3 ! 3 ' 2 GM^ (r) (75) (76) prrip 3(1-/3) One may assume that the contribution of the gas density and mass to the geometric density and geometric mass, respectively, can be neglected. The latter approximations are very well supported by astrophysical observations, which show that the gas represents only a small fraction of the total mass [70ll72]. Therefore, we obtain «A r " 2 , and GM.ir)*^* AnGp^r) « ^JllLJL r ~\ and GM<t> (r) a ^JULJL r> (77) respectively. One may also estimate an upper bound for the cutoff of the geometric mass. The idea is to consider the point at which the decaying density profile of the geometric density associated to the galaxy cluster becomes smaller than the average energy density of the Universe. Let the value of the coordinate radius at the point where the two densities are equal to be nf r . Then at this point p,p{R? r ) = p un iv, where p un iv is the mean energy density of the universe. By assuming p un iv = Pc = 3H 2 /8irG — 4.6975 x 10 _30 /i| g/cm~ 3 , where H = 5O/150 km/Mpc/s J7CJ, we obtain The total geometric mass corresponding to this value is M™ = M, « } ) = 4.83 x 10 16 /5 3 / 2 (|gf ^ K$M®. (79) This value of the mass is consistent with the observations of the mass distribution in the clusters of galaxies. However, according to f(R) modified theories of gravity, we predict that the geometric mass and its effects extends beyond the virial radius of the clusters, which is of the order of only a few Mpc. By assuming that R^ ~ -Rl , we obtain the following relation between the virial and the baryonic mass of the cluster For a cluster with gas temperature T g — 5x 10 7 K, (3 = 1/2 and Ry = 2 Mpc we obtain My ~ 32M, a relation which is consistent with the astronomical observations 1701 . 14 IV. SUMMARY AND DISCUSSION Cosmology has entered a 'golden age', in which the rapid development of increasingly high-precision data has turned it from a speculative to an observationally based science. Recent experiments call upon state of the art technology to provide detailed information about the contents and history of the Universe. These experiments include the Hubble Space Telescope, the NASA WMAP satellite instrument, that measures the temperature and polarisation of the CMBR, and the Sloan Digital Sky Survey (SDSS), that is automatically mapping the properties and distribution of 1 million galaxies. High-precision cosmology has allowed us to tie down the parameters that describe our Universe with growing accuracy. The standard model of cosmology is remarkably successful in accounting for the observed features of the Universe. However, there remain a number of fundamental open questions at the foundation of the standard model. In particular, we lack a fundamental understanding of the acceleration of the late universe. Recent observations of supernovae, together with the WMAP and SDSS data, lead to the remarkable conclusion that our universe is not just expanding, but has begun to accelerate [l|. What is the so-called 'dark energy' that is driving the acceleration of the universe? Is it a vacuum energy or a dynamical held ("quintessence")? Or is the acceleration due to infra-red modifications of Einstein's theory of General Relativity? How is structure formation affected in these alternative scenarios? What will the outcome be of this acceleration for the future fate of the universe? The aspects of these fundamental questions whose resolution is so important for theoretical cosmology, need to look beyond the standard theory of gravity. It is clear that these questions involve not only gravity, but also particle physics. String theory provides a synthesis of these two parts of physics and is widely believed to be moving towards a viable quantum gravity theory. One of the key predictions of string theory is the existence of extra spatial dimensions. In the brane-world scenario, motivated by recent developments in string theory, the observed 3-dimensional universe is embedded in a higher-dimensional spacetime [231 . The generalization of the Dvali-Gabadadze-Porrati (DGP) brane models [3l| lead to late-accelerating cosmologies |32|], even in the absence of a dark energy field. This exciting feature of "self acceleration" may help towards a new resolution to the dark energy problem, although this model deserves further investigation as a viable cosmological model [73|]. It will be interesting to generalize the DGP model with the inclusion of a Gauss-Bonnet (GB) [33[, and it will also be important to investigate the effects of the GB term relatively to the issues of strong coupling and ghosts in the DGP models. Infra-red modifications to General Relativity, where the consistency of various candidate models, including 4-dimensional modifications to the Einstein-Hilbert action, especially GB modifications with a scalar field coupling, have also be analyzed [57] . In this context, a more general modification of the Einstein-Hilbert gravitational Lagrangian density in the form of L — f(R) has recently been extensively analyzed. Relatively to the construction of "quintessential" dark energy models, recent fits to observational data indicate that an evolving equation of state crossing the phantom divide is mildly favored [7J, UM- ^ n a cosmological setting, it has also been shown that the transition into the phantom regime, a mixture of various interacting non-ideal fluids is necessary [6] , with important implications to the model construction of dark energy. If confirmed in the future, this behaviour holds important implications to the model construction of dark energy. The latter models, considering a redshift dependent equation of state, possibly provide better fits to the most recent and reliable SN la supernovae Gold dataset. Deciding between these possible sources of the cosmic acceleration will be one of the major objectives in cosmology in the next decade with several surveys and experiments to address the nature of dark energy. One may mention new several major SNIa supernovae projects, such as the SuperNova Legacy Survey (SNLS), SDSS-II, Destiny, the Large Synoptic Survey Telescope (LSST), Dark Energy Survey (DES) and the SuperNova Acceleration Probe (SNAP). Other dark energy probes include the Dark UNiverse Explorer (DUNE); the Wide Field Multi-Object Spectrograph (WFMOS), which will perform surveys to measure dark energy and the history of our Galaxy; and the Panoramic Survey Telescope and Rapid Response System (PanSTARRS), amongst others. All of these aspects present an extremely fascinating aspect for the above-mentioned experiments and for future theoretical research. Acknowledgments I thank Orfeu Bertolami, Christian Boehmer, Tiberiu Harko, Kazuya Koyama, Roy Maartens and Antonios Papa- zoglou for helpful comments, and acknowledge funding by Fundagao para a Ciencia e a Tecnologia (FCT)-Portugal through the grant SFRH/BPD/26269/2006. 15 [1] Perlmutter, S. et al. 1999, Astrophys. J., 517, 565; Riess, A. G. et al. 1998, Astron. J., 116, 1009; Riess, A. G. et al. 2004, Astrophys. J. 607, 665; Grant, A. et al 2001, Astrophys. J., 560, 49; Perlmutter, S., Turner, M.S., and White, M. 1999, Phys. Rev. Lett., 83, 670; Bennett, C.L. et al 2003, Astrophys. J. Sup pl., 148, 1; Hinshaw, G. et al, arXiv:astro-ph/0302217 . [2] Copeland, E.J., Sami, M., and Tsujikawa, S. 2006, Int. J. Mod. Phys. D, 15, 1753. [3] Douspis, M., Riazuelo, A., Zolnierowski, Y., and Blanchard, A. 2003, Astrophys. J. Suppl., 148, 135; Wang, L., Caldwell , R.R., Ostriker, J .P., and Steinhardt, P.J. 2000, Astrophys. J., 530, 17. [4] Turner, M.S. 2001, |astro-ph /0108103 [5] Caldwell, R.R. 2002, Phys. Lett. B, 545, 23; Caldwell, R.R., Kamionkowski, M., and Weinberg, N.N. 2003, Phys. Rev. Lett., 91, 071301. [6] Vikman, A. 2005, Phys. Rev. D, 71, 023515. [7] Guo, Z., Piao, Y., Zhang, X., and Zhang, Y. 2005, Phys. Lett. B, 608, 177; Zhang, X.F., Li, H., Piao, Y.S., and Zhang, X.M. 2006, Mod. Phys. Lett. A, 21, 231; Perivolaropoulos, L. 2005, Phys.Rev. D, 71, 063503; Wei, FL, and Cai, R.G. 2005, Class. Quant. Grav., 22, 3189; Li, M.z., Feng, B., and Zhang, X.m. 2005, J. Cosmol. Astropart. Phys., 0512, 002; Stefancic, H. 2005, Phys. Rev. D, 71, 124036; Anisimov, A., Babichev, E., and Vikman, A. 2005, J. Cosmol. Astropart. Phys., 0506, 006; Wang, B., Gong, Y.g., and Abdalla, E. 2005, Phys. Lett. B, 624, 141; Nojiri, S., and Odintsov, S.D. 2006, Gen. Rel. Grav., 38, 1285; Aref'eva, I.Y., Koshelev, A.S., and Vernov, S.Y. 2005, Phys. Rev. D, 72, 064017; Zhao, G.B., Xia, J.Q., Li, M., Feng, B., and Zhang, X. 2005, Phys. Rev. D, 72, 123515; Tsujikawa, S. 2005, Phys. Rev. D, 72, 083512. [8] Kamenshchik, AY., Moschella, U., and Pasquier, V. 2001, Phys. Lett. B 511, 265; Bento, M.C., Bertolami, O., and Sen, A.A. 2002, Phys. Rev. D 66, 043507; Bento, M.C., Bertolami, O., and Sen, A.A. 2003, Phys. Rev. D 67, 063003; Bento, M.C., Bertolami, O., and Sen, A.A. 2003, Gen. Rel. Grav. 35, 2063; Bento, M.C., Bertolami, O., and Sen, A.A. 2003, Phys. Lett. B 575, 172; Amendola, L., Finelli, F., Burigana, C, and Carturan, D. 2003, JCAP 0307, 005; Bento, M.C., Bertolami, O., and Sen, A.A. 2004, Phys. Rev. D 70, 083519. [9] Bilic, N., Tupper, G.B., and Viollier, R.D. 2002, Phys. Lett. B 535, 17. [10] Weyl, H. 1921, Space, Time, Matter, Chapter IV, New York, Dover; Eddington, A. 1924, The Mathematical Theory of Relativity, Chapter IV, London, CUP; Lanczos, K. 1938, Ann. Math., 39, 842; Buchdahl, H.A. 1948, Proc. Edin. Math. Soc, 8, 89; Pais, A., and Uhlenbeck, G.E. 1950, Phys. Rev., 79, 145; Utiyama, R., and de Witt, B. 1962, J. Math. Phys., 3, 608; Havas, P. 1977, Gen. Rel. Grav., 8, 631; Stelle, K. 1978, Gen. Rel. Grav., 5, 353; Mannheim, P.D., and Kazanas, D. 1989, Astrophys. J., 342, 635; Kazanas, D., and Mannheim, P.D. 1991, Astrophys. J. Suppl., 76, 421. [11] Buchdahl, H.A. 1970, Mon. Not. Roy. Astron. Soc, 150, 1. [12] Kerner, R. 1982, Gen. Rel. Grav., 14, 453. [13] Duruisseau, J. P., Kerner, R., and Eysseric, P. 1983, Gen. Rel. Grav., 15, 797. [14] Barrow, J.D., and Ottewill, A.C. 1983, J. Phys. A: Math. Gen., 16, 2757. [15] Starobinsky, A.A. 1980, Phys. Lett. B, 91, 99. [16] Carroll, S.M., Duvvuri, V., Trodden, M., and Turner, M.S. 2004, Phys. Rev. D, 70, 043528. [17] Amendola, L., Polarski, D., and Tsujikawa, S. 2007, Phys. Rev. Lett., 98, 131302; Capozziello, S., Nojiri, S., Odintsov, S.D., and Troisi, A. 2006, Phys. Lett., B, 639, 135; Nojiri, S., and Odintsov, S.D. 2006, Phys. Rev. D, 74, 086005; Amarzguioui, M., Elgaroy, O., Mota, D.F., and Multamaki, T. 2006, Astron. Astrophys., 454, 707; Amendola, L., Gannouji, R., Polarski, D., and Tsujikawa, S. 2007, Phys. Rev. D, 75, 083504; Koivisto, T. 2007, Phys. Rev. D, 76, 043527; Starobinsky, A.A. 2007, JETP Lett., 86, 157; Li, B., Barrow, J.D., and Mota, D.F. 2007, Phys. Rev. D, 76, 044027; Perez Bergliaffa, S.E. 2006, Phys. Lett. B, 642, 311; Santos, J., Alcaniz, J.S., Reboucas, M.J., and Carvalho, F.C. 2007, Phys. Rev. D, 76, 083513; Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D., and Zerbini, S. 2005, J. Cosmol. Astropart. Phys., 0502, 010; 16 Faraoni, V. 2005, Phys. Rev. D, 72, 061501; Faraoni, V. 2005, Phys Rev. D, 72 , 124005; Sokolowski, L.M. 2007, |arXiv:gr-qc/0702097| Cognola, G., Gastaldi, M., and Zerbini, S. 2007, |arXiv:gr-qc/0701138| Bohmer, C.G., Hollenstein, L., and Lobo, F.S.N. 2007, Phys. Rev. D, 76, 084005; Carloni, S., Dunsby, P.K.S., and Troisi, A. 2007, arXiv:gr-qc/0707.0106; Ananda, K.N., Carloni, S., and Dunsby, P.K.S. 2007, arXiv:gr-qc/0708.2258; Capozziello, S., Cianci, R., Stornaiolo, C, and Vignolo, S. 2007, Class. Quant. Grav., 24, 6417; Tsujikawa, S. 2008, Phys. Rev. D, 77, 023507; Nojiri, S., Odintsov, S.D., and Tretyakov, PV. 2007, Phys. Lett. B, 651, 224; Nojiri, S., and Odintsov, S.D. 2007, Phys. Lett. B, 652, 343; Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D., Sebastiani, L., and Zerbini, S. 2008, Phys. Rev. D, 77, 046009. [18] Chiba, T. 2003, Phys. Lett. B, 575, 1; Erickcek, A.L., Smith, T.L., and Kamionkowski, M. 2006, Phys. Rev. D, 74, 121501; Chiba, T., Smith, T.L., and Erickcek, A.L. 2007, Phys. Rev. D, 75, 124014; Nojiri, S., and Odintsov, S.D. 2008, Phys. Lett. B, 659, 821; Capozziello, S., Stabile, A., and Troisi, A. 2007, Phys. Rev. D, 76, 104019; Capozziello, S., Stabile, A., and Troisi, A. 2008, Class. Quant. Grav., 25, 085004. [19] Olmo, G.J. 2007, Phys. Rev. D, 75, 023511. [20] Hu, W., and Sawicki, I. 2007, Phys. Rev. D, 76, 064004; Nojiri, S., and Odintsov, S.D. 2003, Phys. Rev. D, 68, 123512; Faraoni, V. 2006, Phys. Rev. D, 74, 023529; Faulkner, T., Tegmark, M., Bunn, E.F., and Mao, Y. 2007, Phys. Rev. D, 76, 063505; Zhang, P.J. 2007, Phys. Rev. D, 76, 024007; Capozziello, S., and Tsujikawa, S. 2008, Phys. Rev. D, 77, 107501; Sawicki, I., and Hu, W. 2007, Phys. Rev. D, 75, 127502; Amendola, L., and Tsujikawa, S. 2008, Phys. Lett. B, 660, 125. [21] Capozziello, S., Cardone, V.F., and Troisi, A. 2006, J. Cosmol. Astropart. Phys., 0608, 001; Capozziello, S., Cardone, V. F., and Troisi, A. 2007, Mon. Not. R. Astron. Soc, 375, 1423. [22] Borowiec, A., Godlowski, W., and Szydlowski, M. 2007, Int. J. Geom. Meth. Mod. Phys., 4, 183 [23] Martins, C.F., and Salucci, P. 2007, Mon. Not. Roy. Astron. Soc, 381, 1103. [24] Boehmer, C.G., Harko, T., and Lobo, F.S.N. 2008, arXiv:gr-qc/0709.0046. [25] Boehmer, C.G., Harko, T., and Lobo, F.S.N. 2008, J. Cosmol. Astropart. Phys., 0803, 024. [26] Bertolami, O., Boehmer, C.G., Harko, T., and Lobo, F.S.N. 2007, Phys. Rev. D, 75, 104016. [27] Bertolami, O., and Paramos, J. 2008, Phys. Rev. D, 77, 084018. [28] Binney, J., and Tremaine, S. 1987, Galactic Dynamics, Princeton, Princeton University Press; Persic, M., Salucci, P., and Stel, F. 1996, Month. Not. R. Astron. Soc, 281, 27; Borriello, A., and Salucci, P. 2001, Month. Not. R. Astron. Soc, 323, 285; Salucci, P., Lapi, A., Tonini, C, Gentile, G., Yegorova, I., and Klein, U. 2007, Mon. Not. Roy. Astron. Soc, 378, 41. [29] Maartens, R. 2004, Living Rev. Rel. 7, 7. [30] Randall, L., and Sundrum, R. 1999, Phys. Rev. Lett., 83, 3370; Randall, L., and Sundrum, R. 1999, Phys. Rev. Lett., 83, 4690. [31] Dvali, G., Gabadadze, G., and Porrati, M. 2000, Phys. Lett. B, 485, 208. [32] Deffayet, C. 2001, Phys. Lett. B, 502, 199. [33] Brown, R.A., Maartens, R., Papantonopoulos, E., and Zamarias, V. 2005, J. Cosmol. Astropart. Phys., 11, 008. [34] Lue, A. 2006, Phys. Rep., 423, 1. [35] Koyama, K., and Maartens, R. 2006, J. Cosmol. Astropart. Phys., 0601, 016. [36] Lue, A. 2003, Phys. Rev. D, 67, 064004. [37] Burgess, C.P. 2004, Ann. Phys. 313, 283. [38] Sotiriou, T.P., and Faraoni, V. 2008, arXiv:gr-qc/0805.1726. [39] Nojiri, S., and Odintsov, S.D 2007. Int. J. Geom. Meth. Mod. Phys., 4, 115; [40] Koyama, K. 2008, Gen. Rel. Grav., 40, 421. [41] Buchdahl, H.A. 1970, Month. Not. R. Astron. Soc, 150 , 1; Barrow, J.D., and Ottewill, A.C. 1983, J. Phys. A: Math. Gen., 16, 2757. [42] Ferraris, M., Francaviglia, M., and Volovich, I. 1993, |arXivlgr^qc /9303007 Vollick, D.N. 2003, Phys. Rev. D, 68, 063510; Flanagan, E.E. 2003, Class. Quant. Grav., 21, 417; Meng, X.H., and Wang, P. 2004, Phys. Lett. B, 584, 1; Li, B., and Chu, M.C. 2006, Phys. Rev. D, 74, 104010; Poplawski, N.J. 2006, Phys. Rev. D, 74, 084032; Li, B., Chan, K.C., and Chu, M.C. 2007, Phys. Rev. D, 76, 024002; Li, B., Barrow, J.D., and Mota, D.F. 2007, Phys. Rev. D, 76, 104047; Iglesias, A., Kaloper, N., Padilla, A., and Park, M. 2007, Phys. Rev. D, 76, 104001. [43] Sotiriou, T.P., and Liberati, S. 2007, Ann. Phys., 322, 935. 17 [44] Bertolami, O., and Paramos, J. 2008, arXiv:gr-qc/0805.1241. [45] Faraoni, V. 2007, Phys. Rev. D, 76, 127501. [46] Nojiri, S., and Odintsov, S.D. 2004, Phys. Lett. B, 599, 137; Nojiri, S., and Odintsov, S.D. 2004, Proceedings of Science, WC2004, 024; Allemandi, G., Borowiec, A., Francaviglia, M., and Odintsov, S.D. 2005, Phys. Rev. D, 72, 063505. [47] Mukohyama, S., and Randall, L. 2004, Phys. Rev. Lett., 92, 211302. [48] Koivisto, T. 2006, Class. Quant. Grav., 23, 4289. [49] Sotiriou, T.P., and Faraoni, V. 2008, arXiv:gr-qc/0805.1249. [50] Schutz, B.F. 1970, Phys. Rev. D, 2, 2762. [51] Brown, J.D. 1993, Class. Quant. Grav., 10, 1579. [52] Bertolami, O., Lobo, F.S.N., and Paramos, J. 2008, arXiv:gr-qc/0806.4434. [53] Hawking, S.W., and Ellis, G.F.R. 1973, The Large Scale Structure of Spacetime, (Cambridge University Press, Cambridge). [54] Capozziello, S., Carloni, S., and Troisi, A. 2003, Recent Res. Dev. Astron. Astrophys., 1, 625. [55] Sotiriou, T.P. 2007, arXiv:gr-qc/0710.4438. [56] Nojiri, S., and Odintsov, S.D. 2004, Gen. Rel. Grav., 36, 1765. [57] Nojiri, S., Odintsov, S.D., and Sasaki, M. 2005, Phys. Rev. D, 71, 123509. [58] Koivisto, T., and Mota, D.F. 2007, Phys. Rev. D 75, 023518; Koivisto, T., and Mota, D.F. 2007, Phys. Lett. B 644, 104. [59] Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D., and Zerbini, S. 2006, Phys. Rev. D, 73, 084007; Nojiri, S., Odintsov, S.D., and Gorbunova, O.G. 2006, J. Phys. A: Math. Gen., 39, 6627. [60] Nojiri, S., and Odintsov, S.D. 2005, Phys. Lett. B, 631, 1. [61] Koyama, K. 2007, Class. Quant. Grav., 24, R231. [62] Koyama, K., and Maartens, R. 2006, J. Cosmol. Astropart. Phys., 0601, 016. [63] Maartens, R., and Majerotto, E. 2006, Phys. Rev. D 74, 023004; Bento, M.C., Bertolami, O., Reboucas, M.J., and Santos, N.M.C. 2006, Phys. Rev. D 73, 103521. [64] Koyama, K. 2005, Phys. Rev. D, 72, 123511; Gorbunov, D., Koyama, K., and Sibiryakov, S. 2006, Phys. Rev. D, 73, 044016; Izumi, K., Koyama, K., and Tanaka, T. 2007, J. High Energy Phys., 0704, 053. [65] de Rham, C, Dvali, G., Hofmann, S., Khoury, J., Pujolas, O., Redi, M., and Tolley, A.J. 2007, arXiv:hep-th/0711.2072. [66] Overduin, J.M., and Wesson, P.S. 2004, Phys. Rep., 402, 267; Bohmer C.G., and Harko, T. 2007, Month. Not. R. Astron. Soc, 379, 393; Sin, S.J. 1994, Phys. Rev. D, 50, 3650; Silverman, M.P., and Mallett, R.L. 2002, Gen. Rel. Grav., 34, 633; Ferrer, F., and Grifols, J.A. 2004, J. Cosmol. Astropart. Phys., 0412, 012; Bohmer C.G., and Harko, T. 2007, J. Cosmol. Astropart. Phys., 06, 025. [67] Milgrom, M. 1983, Astrophys. J., 270, 365; Sanders, R.H. 1984, Astron. Astrophys., 136, L21; Bekenstein, J.D. 2004, Phys. Rev. D, 70, 083509; Moffat, J.W., and Sokolov, I.Y. 1996, Phys. Lett. B, 378, 59; Mannheim, P.D. 1997, Astrophys. J., 479, 659; Brownstein, J.R., and Moffat, J.W., Roberts, M.D. 2004, Gen. Rel. Grav., 36, 2423; Brownstein, J.R., and Moffat, J.W. 2006, Astrophys. J., 636, 721; Brownstein, J.R., and Moffat, J.W. 2006, Month. Not. R. Astron. Soc, 367, 527; Mak, M.K., and Harko, T. 2004, Phys. Rev. D, 70, 024010; Harko, T., and Cheng, K.S. 2006, Astrophys. J., 636, 8; Bohmer, C.G., and Harko, T. 2007, Class. Quantum Grav., 24, 3191; Harko, T., and Cheng, K.S. 2007, Phys. Rev. D, 76, 044013. [68] Lindquist, R.W. 1966, Ann. Phys, 37, 487; Maartens, R., and Maharaj, S.D. 1985, J. Math. Phys., 26, 2869; Bildhauer, S. 1989, Class. Quantum Grav., 6, 1171; Banach, Z., and Piekarski, S. 1994, J. Math. Phys., 35, 4809. [69] Jackson, J.C. 1970, Month. Not. R. Astr. Soc, 148, 249. [70] Reiprich, T.H., and Boringer, H. 2002, Astrophys. J., 567, 716. [71] Harko, T., and Cheng, K.S. 2007, Phys. Rev. D, 76, 044013. [72] Arnaud, M. 2005, "X-ray observations of clusters of galaxies", in "Background Microwave Radiation and Intracluster Cosmology", Proceedings of the International School of Physics "Enrico Fermi", edited by F. Melchiorri and Y. Rephaeli, published by IOS Press, The Netherlands, and Societa Italiana di Fisica, Bologna, Italy, p. 77; Schuecker, P., Bohringer, H., Arzner, K., and Reiprich, T.H. 2001, Astron. Astrophys., 370, 715; Carlberg, R.G., Yee, H.K.C., and Ellingson, E. 1997, Astrophys. J., 478, 462. [73] Kofinas, C, Maartens, R., and Papantonopoulos, E. 2002, J. High Energy Phys., 0310, 066. [74] Feng, B., Wang, X., and Zhang, X. 2005, Phys. Lett. B, 607, 35. [75] Upadhye, A., Ishak, M., and Steinhardt, P.J. 2005, Phys. Rev. D, 72, 063501.